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Abstract 

Using a sample of almost 7000 strong Mgll absorbers with 0.4 < z < 2.2 detected in the SDSS 
DR4 dataset, we investigate the gravitational lcnsing and dust extinction effects they induce on 
background quasars. After carefully quantifying several selection biases, we isolate the reddening 
effects as a function of redshift and absorber rest equivalent width, Wo. We find the amount of dust 
to increase with cosmic time as r(z) oc (1 + z) _11±0 ' 4 , following the evolution of cosmic star density 
or integrated star formation rate. We measure the reddening effects over a factor 30 in E(B-V) and 
we find that r oc (Wo) 1 , providing us with an important scaling for theoretical modeling of metal 
absorbers. We also measure the dust-to-metals ratio and find it similar to that of the Milky Way. 
In contrast to previous studies, we do not detect any gravitational magnification by Mgll systems. 
We measure the upper limit /x < 1.10 and discuss the origin of the discrepancy. Finally, we estimate 
the fraction of absorbers missed due to extinction effects and show that it rises from 1 to 50% in the 
range 1 < Wq < 6 A. We parametrize this effect and provide a correction for recovering the intrinsic 
dN/dWo distribution. 

Subject headings: quasars - absorbers: Mgll - gravitational lensing: statistical - reddening - dust 



1. INTRODUCTION 

Quasar absorption lines provide us with a unique tool 
to probe the gas content in the Universe. They offer 
an unmatched sensitivity up to high redshift, allow us to 
detect arbitrarily faint objects, and constrain the gas dis- 
tribution in velocity space. However, the lack of direct 
spatial information tends to be a limitation for under- 
standing the nature of these systems and fundamental 
questions still remain. One way to overcome this problem 
is to relate the gas distribution to that of other compo- 
nents such as dark matter, stars and dust by measuring 
their cross-correlations. 

Galactic environments have been studied using various 
types of absorption lines. Among metal lines, the Mgll 
doublet, AA2796, 2803 A, has been extensively used due 
to the strength and wavelength of the transitions mak- 
ing i t easily detecta ble from the ground at 0.4 5, z < 
2.2 (lLanzetta et all (fl987h: ISteidel fc Sargentl { 1992t h 
iNestor et all (|2005h ; [Prochter et all (|2006t )). Mgll turns 
out to be a sensitive tracer of galactic environment: it 
arises in gas spanning more than five decades of neu - 
tral hydrogen column density (e.g. IChurchill et al.ll2000f ). 
and it has been shown that the detection of a strong Mgll 
absorber indicates the presence of a nearby ~ 0.1 — 10 L* 
galaxy along the l i ne-of-sight (|Bergeron fc Boiss e|[T991t 
Steide etldlll994t iBergeron et al.l 11992b iLe Brun et all 



1993t ISteidel et al.l [19941. Il997f). and at low redshift 
Rao et al.l (|2006f )~have shown that strong Mgll absorbers 
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are good tracers of Damped Lyman alpha systems. For 
reviews of t he progress on t he de termi nation of these 
system s, see iTripp fc Bowenl (|2005h and IChurchill etldl 

pool . 

Recently, the size and homogeneity the SDSS 
has allowed severa l statistical analyzes to be done: 
iBouche et al.1 ()2006l ) measured the cross-correlation be- 
tween z ~ 0.5 luminous r ed galaxies and Mgll absorbers. 
IZibetti et"aH (|2005l 120071 ) investigated the distribution of 
absorber impact parameters and the emission properties 
of Mgll related galaxies. They have shown that a cross- 
correlation between light and gas can be detected up to 
200 kpc and that stronger absorbers are related to bluer 
galaxies. The SDSS has also provided us with closely 
separated quasar pairs allowing us to analyze the spec- 
trum of background quasars to probe the gas d i stribu - 
tion around foreground ones (|Bowen et al.ll2(M 12001 . 
Such a technique allows us to compare the gaseous halos 
of quasars and galaxies and it can ultimately be used 
to test the quasar un ification scheme on large scales 
(jChelouche et al.ll2007T i. 

However, despite significant observational progress 
fundamental questions such as the nature of the gas seen 
in absorption and its spatial distribution still remain. 
More importantly, we still do not understand the pro- 
cesses) responsible for the observed values of Mgll rest 
equivalent widths. As strong Mgll absorbers are known 
to be tracers of galaxies, they are expected to induce 
gravitational lensing and reddening effects to the light of 
their background quasars. Correlations between quasar 
magnitudes and the presence or strength of metal lines 
can therefore be used in order to quantify the relation- 
ships between gas, dark matter and dust. Such an ap- 
proach also allows us to quantify the extinction and mag- 
nification biases in samples of optically selected quasars. 

1.1. Dust reddening 
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The presence of dust associated with variou s types of 
absorb er s has been reported by many authors: iPei et all 
(|1991h & (|1992ft found quasars with DLAs to be on 
averag e redder than quasars without. Richard s et al.l 
(2001) observed a similar ef fect with CIV abs o rption 
lines. Using the 2dF survey IMenard fc Perouxl (2003) 
found evidence of re d denin g due to strong Mgll sys- 
tems. Hopk ins et al.l (|2004h showed that a fraction of 
SDSS quasars are reddened with SMC-type dust and ar- 
gued that this effect comes predominantly from dust lo- 
cated at the quasar re d shifts. Using the S PSS database, 
iWild fc Hewettl (l27)(^ : IWuTeTall p006l ) fou nd an E(B- 
V)~ .06 associated wi t h Cal l absorbers, and I York et al.l 
( 2006) & iKhare et all (|2005h measured the mean red- 
dening effects of about 800 quasars with z > 1 
Mgll absorbers. Constraints on dust can also be ob- 
tained by measuring the relative abundances of volatile 
and refract ory elements in order t o infer a depletion 
level (e.g. IVladilo ifc Peroul[2fM iVladilo et all [20061 : 
iNestor et alT~ll2003f ). 

While the presence of dust reddening associated with 
metal absorbers has been convincingly shown, its de- 
tailed properties (dependencies on redshift, rest equiv- 
alent width, etc.) need to be quantified in order to un- 
derstand the underlying physics and connect theory to 
observations. 

The presence of dust along given lines-of-sight also im- 
plies the existence of a dust obscuration bias that can 
affect quasar and abs orber number count s. Using 75 
radio-selected quasars, lEllison et alJ (|2004f ) investigated 
such an effect and did not find any difference (at the 
20% level) between the incidence of Mgll absorbers in 
radio-selected and optically-selected quasars. This result 
shows that in general optically-selected quasars are not 
strongly affected by an extinction bias, but it does not 
provide us with the full information: it does not show 
how the shape of the observed distribution of absorber 
rest equivalent widths is affected by dust extinction. 

In order to unveil the relationships between extinc- 
tion and absorption properties, we present a n ew analysis 
based on the Sloan Digital Sky Survey (SDSS. lYork et all 
2000) using the fourth data release (DR4). After having 
quantified the efficiency of the SDSS pipeline to detect 
reddened quasars, we use a sample of about 7,000 Mgll 
absorbers - i.e. almost an order of magnitude larger 
than previous studies - to measure the excess reddening 
induced by these systems and investigate its dependence 
on absorption rest equivalent width and redshift. 

1.2. Gravitational lensing 

In addition to cause reddening effects, the presence of 
a galaxy along the line-of-sight is expected to give rise to 
gravitational lensing effects. Several cases of strong grav- 
itational lensing o f quasars due to an a bsorbing galaxy 
are known (e.g. Turnshck et alJ Il997t ). However, in 
general, the impact parameters of Mgll absorbers are 
larger than a few kpc and gravitational lensing effects 
occur in the weak regime and solely c hange quasa r mag- 
nitudes (IBartelmann fc Loebl 1199a: iPema et al.l Il997t 
Smette et al.1 Il997t iBartelmann fc Loebl 119981: IMenard! 



2005). The detection of such an effect would ultimately 



allow us to constrain the mass of Mgll systems. 

Several authors have attempted to detect the 
statistical magnification induced by Mgll systems. 



Some have looked for redshift distribution changes 
(jThomas fc Websterl 119901 : iBorgeest fc Mehlertl Il993| ) 
but did not find any compelling evidence for gravi- 
tational lensing. Others have investigated the occur- 
rence of metal ab sorbers (CIV, SilV, Mgll) on bright 
and faint quasars (lYork et al.lll991t IVanden Berk et all 
Il996t iRichards et al.lll999f) and reported possible mag- 
nification effects for c ertain absorber species only. 
IMenard fc Perouxl (|2003f ) used the 2dF-Quasar survey to 
compare the number of quasars with and without Mgll 
absorbers as a function of magnitude and reported a 
relative excess of optically bright radio-selected quasars 
with absorbers as well as indications of reddening ef- 
fects. Murphy & Liske (2004) applied the same method 
to a sample of Damped Lyman-a systems and reported a 
similar trend found at 2a. Using radio-selected quasars, 
lEllison et al.l (|2004h reported a possible excess of opti- 
cally bright quasars wi t h Mg ll absorbers, and more re- 
cently iProchaska et al.l (|2005h found that quasars with 
high-redshift DLAs are on average brighter than quasars 
without. The various results found among all these an- 
alyzes turn out to be rather difficult to combine be- 
cause the absorber species, redshift and equivalent width 
ranges differ. In addition, as noted in some of these ana- 
lyzes, biases due to selection effects can mimic the lensing 
signal: it is easier to detect absorption lines in high qual- 
ity (or high S/N) spectra, however, these spectra tend to 
correspond to bright objects. This selection bias results 
in a preference for detecting absorption lines in brighter 
quasars, which is an effect similar to the signal of inter- 
est. Unfortunately, none of these works could convinc- 
ingly demonstrate that the level of such systematics was 
sufficiently low in their analysis, leaving the amplitude 
of lensing effects rather uncertain. 

In this paper we attempt to make a significant improve- 
ment in making such a measurement by overcoming the 
problem of systematic effects using Monte Carlo simula- 
tions and by reaching a much higher precision thanks to 
the size and homogeneity of the SDSS. We use a robust 
semi-automatic absorption line finder algorithm, care- 
fully select a population of reference quasars, test poten- 
tial biases, and use Monte Carlo simulations to measure 
and quantify the level of systematics of our selection pro- 
cedure. To detect a magnification si gnal we will apply the 
method proposed by IMenard] (|2005f ) and look for a corre- 
lation between the presence of absorbers and a change in 
the mean magnitude of their background quasars. Con- 
trarily to previous studies, this approach allows the use 
of a non-parametric estimator. 

The outline of the paper is as follows: in $2] we intro- 
duce the formalism of dust extinction and gravitational 
magnification. We describe the data used in our analysis 
in |J3] as well as the selection of the different samples of 
quasars. The magnitude shifts and color changes due to 
the presence of absorbers along quasar lines-of-sight are 
presented in ||4] We discuss the results and conclude in 

m 

2. FORMALISM 

Strong Mgll absorbers (usually defined with 
Wo (2796) > 0.3 A) are known to be tracers of galaxies 
with luminosities ranging from 0.1 to several L*. The 
presence of such a system seen along the line-of-sight of 
a background source can therefore modify its brightness 
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in two ways: first, the presence of dust in and/or around 
the absorber can extinct and redden the source's light 
and second, it can act as a gravitational lens and amplify 
the flux of the source. This can be described by 

f(A)=f ref (A)Ate-^, (1) 

where f(A) is the observed flux of a source behind an ab- 
sorber, f re f (A) is the flux that would be observed without 
intervening system, t\ is the optical depth for extinction 
induced by the galaxy and fi the gravitational amplifica- 
tion. The corresponding magnitude change can then be 
written 

dm = m — m rc f 

= -2.5 1og(A0 + j^7*. (2) 

These two contributions have different wavelength and 
rcdshift contributions. One can therefore attempt to con- 
strain them separately. 

2.1. Measuring reddening effects 

Following [Menardl (|2005| ), we introduce the observed 
mean magnitude difference between a population of 
quasars with absorbers and one without: 

A(mj) = (mj) - (m^ re f) , (3) 

where nij is the magnitude in a given band j. In the 
next section we will measure this quantity using a large 
sample of Mgll absorbers from the SDSS. We will then 
investigate how A(mj) varies as a function of absorber 
rest equivalent width and redshift. Similarly, we can then 
define the mean excess color 

(E(J — k)) = A(m J ) — A(mfc) (4) 

where j and k denote two different bands. It should 
be noted that such a quantity is not sensitive to grav- 
itational lensing effects as it no longer depends on the 
magnification fi. It can be used to probe the statistical 
properties of the extinction curve of absorbers systems 
as well as the dust column densities as a function of 
absorber rest equivalent width. 

In optical surveys, sources like quasars are both mag- 
nitude and color selected. When dust extinction occurs, 
such sources can become undetectable because they be- 
come too faint and/or too red. As long as sources are no 
longer detectable because of a limiting magnitude and 
not because of a color cut, we have 

(E(j-k)) = ^{r(X j )-r{X k )). (5) 

However, if the sources become undetectable because 
they no longer satisfy a color cut, the measured redden- 
ing is then lower than the intrinsic one. Such a property 
depends on the details of the quasar selection procedure 
and will be addressed in section 13.21 

The current study is based on the photometry of SDSS 
quasars using the u,g,r,i,z filter set. However, his- 
torically people have quantified color excesses using the 
Johnson B and V filters. In order to be able to present 
our results in similar units, we provide here a list of use- 
ful conversions. First, it should be noted that previ- 
ous statistical studies (Kh are et a.1.1 l2005t iMenard et al.l 



l2005t lYork etldl 12001 iWild et al.ll200l have indicated 
that, on average, the extinction curves of Mgll selected 
galaxies are consistent with that of the SMC, i.e. do 
not present an excess around the rest frame wavelength 
of 0.2^m (in the optical, such a feature can only be de- 
tected for systems at redshift greater than ~ 1 and is 
therefore irrelevant for lower redshift galaxies). In this 
paper we will therefore consider only SMC-type extinc- 
tion curves. We will also see in section $4] that such a 
property is confirmed by the present analysis. 

In the range 0.1 to 1 fim, the SMC extinction curve 
6 can be very well appro ximated by a power law Ay oc 
A" 12 dPrevot et aLlllQSl . By using the g and i filters, 
the conversion between the two color systems therefore 
reads 

^-1.2 _ ^-1.2 

E (9 ~i)= * 1 2 *-i 2 >< E ( B ~ V) 

~ 1.55 x E(B - V) (6) 

and the average wavelengths for a flat spectrum are given 
by {A g , Ai, A B , Ay} = {468, 748, 435, 555} nm. Using the 
effective wavelength for realistic quasar spectra does not 
change the results significantly. In practice, two frames 
are interesting for considering colors excesses: in order to 
quantify the fraction of sources missed because of an ex- 
tinction bias, observed E(B-V) values are relevant. How- 
ever, for estimating the amount of dust associated with 
an absorber system the rest frame E(B-V) is the quantity 
of interest. For an SMC-like extinction curve we have 

E(B - V) obs = E(B - V) rest (1 + z) 1 - 2 . (7) 

In a given band, we will therefore observe almost four 
times more extinction by an absorber at z = 2 than the 
same system located at z — 0. For an SMC extinction 
curve we also have A v ~ Ca x E(B — V) where the pro- 
portionality coefficient Ca depends on the redshift of the 
dust and varies from 3.2 to 2.0 when z is in the range 0.5 
to 2.0, with a mean value of 2.6. 

2.2. Measuring gravitational lensing effects 

In this section we show how magnification effects can 
(or cannot) affect the observed magnitude of distant 
sources like quasars. For more details we refer the rea der 
to the original discussion presented in iMenardl (|2005l ). 

Here we will first assume that extinction effects can be 
neglected (which is expected at sufficiently large wave- 
lengths). Let us consider an area of the sky which is large 
enough for the mean magnification to be close to unity. 
In this area, let us consider a population of sources, with 
a fraction / of them lying behind an absorber system. We 
will write the magnitude distribution of the sources be- 
hind absorbers N(m), and N rc f (m) for the non-absorbed 
sources. We then have the following relation between the 
magnitude number counts: 

N(m) = / X J N rcf (m - 8m) P(Sm) d(Sm) , (8) 

where 8m is the induced magnitude shift and P(i5m) is 
the probability density to have a given value of 8m along 
a quasar line-of-sight. Note that in the case of absorber 

6 Templates of the reddening curves can be obtained from 
| htt p : / / www . astro, princet on .edu/ ~draine| 
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systems, / is not known. Observational estimates of the 
number density of absorbers are necessarily affected by 
lensing and/or extinction effects due to the need of back- 
ground quasars to observe absorber systems. Therefore 
the only usable information from the previous equation 
comes from the shape of the distribution N(m), i.e. its 
moments. In this paper we will focus only on the first 
moment, i.e. we will attempt to constrain the mean mag- 
nification and extinction effects by measuring the mean 
magnitude shift induced by the presence of absorber sys- 
tems. 

By defining (m) and (m) re f to be the observed mean 
magnitudes of the two quasar populations, we see that 
the observable mean magnitude shift (Am b s ) (cq. EJ) can 
be used to probe the unknown distribution P(<5m). It is 
important to note that the o&serued-magnitude changes 
(Am b s ) depend on both the magnification effects and 
the slope of the source luminosity distribution: 

(ATO obs ) = F[N(m)] x dm , (9) 

and in general we h ave (Am^hs) < (dm). For SDSS 
quasars with i < 20, iMenardl (|2005[ ) has shown that we 
have (Am bs) — 0.30 dm where 5m is an achromatic mag- 
nitude change. For a limiting magnitude of i — 19.1, we 
find (Am bs) — 0.25 5m. 

It should be noted that if the number of sources as a 
function of luminosity follows a power law then F = 
and magnitude changes cannot be observed whatever the 
value of the induced magnitude shift dm. This absence 
of observable lensing effects can be understood in the 
following way: gravitational magnification makes each 
source brighter but allows also new sources to be de- 
tectable as they become brighter than the limiting mag- 
nitude of a given survey. This latter effect increases the 
number of faint sources and tends to diminish the mean 
flux of the detected objects. For a power law luminos- 
ity distribution, the flux increase due to magnification is 
exactly canceled by the additional faint sources that en- 
ter the sample due to magnification, and lensing effects 
due to absorbers cannot be observed. Therefore, it is 
possible to detect changes in the magnitude distribution 
of a population lensed by absorbers only if its unlensed 
number counts depart from a power law as a function of 
luminosity. Extinction effects will be affected in a similar 
way but the ability to observe sources in different bands 
allows us to measure reddening effects and then infer the 
related extinction. 

3. DATA & SELECTION BIASES 
3.1. The quasar sample 

In order to isolate the magnitude changes due to the 
presence of intervening absorber systems we need to com- 
pare the magnitude of quasars with absorption lines to 
that of a reference quasar population. Defining such a 
reference population without introducing any selection 
bias is a difficult task. Indeed, the detection of absorp- 
tion lines depends on quasar redshift, spectrum S/N, and 
the extinction and reddening effects they might induce 
on their background source. All these effects have to be 
taken into account in order to carefully select a popula- 
tion of reference quasars. 

The first requirement is to work with a well defined 
sample of quasars. Among the large sample of quasars 



available from the SDSS DR4 sample (N -50,000), we 
will use the primary spectroscopic quasar sample, i.e. 
quasars with i < 19.1 (N —29,000), as the selection of 
these objects is based on a well defined series of color 
cuts (Schneide r et al.l [2002) . We will not make use of 
fainter quasars (observed with SDSS serendipity fibers) 
as they are not selected from well defined selection crite- 
ria. Before characterizing the absorber systems used in 
this study, the first step of our analysis will be to investi- 
gate the sensitivity of the SDSS quasar target algorithm 
to reddening and extinction. 

3.2. Detection of reddened quasars 

Evidence for the presence of dust associated with Mgll 
absorbers has already been shown by several authors by 
isolating the reddening and extinction effects they in- 
duce. Since the selectio n of quasars depends on quasar 
magnitudes and colors ([Schneider et al.ll2002l ). such an 
effect might introduce a bias in the object selection and 
it is therefore important to quantify it. In this section we 
now invest igate the sensitivity o f the SDSS quasar target 
algorithm (Richards et al. 20021) to redde ning and extinc- 
tion (similarly to the analysis done by iRichards et al.1 
2003). In order to do so we use Montecarlo-like simu- 
lations to test the ability of the SDSS quasar algorithm 
to detect quasars for different amounts of reddening and 
extinction. We proceed as follows: 

1. We use a sample of quasars without strong Mgll 
absorber (as described in section |3~4|) . consider the 
SMC extinction curve and apply various amounts 
of extinction. We quantify the amount of reddening 
with the E(B — V) excess. 

2. The applied extinction will cause some objects to 
drop out of the input quasar sample as they be- 
come too faint and/or their colors no longer satisfy 
the series of selection criteria of the SDSS quasar 
selection algorithm, considering only non-radio and 
point sources. We reject the objects no longer la- 
beled as quasars by the SDSS pipeline. 

3. For each value of E(B — V) applied, we mea- 
sure the fraction of quasars missed due to limiting 
magnitude and color cuts effects, wc measure the 
mean u, g, r, i and z magnitudes of the re-observed 
quasars and compare them to those of the initial 
sample of unreddened quasars. 

We illustrate some of the results of this procedure in 
Fig. [TJ From top to bottom we show the ability of 
the SDSS quasar selection algorithm to select quasars 
with increasing reddening and show the effects on quasar 
magnitudes. In each case we present the u, g, r, i and ri- 
band magnitude distributions of the original (solid line 
histograms) and reddened samples (filled histograms). 
These panels show that the fraction of quasars missed 
due to dust extinction becomes substantial for observed 
reddening values of E(B — V) larger than 0.2. In this 
figure we also show the mean magnitude of each sample 
with a vertical line. As the amount of extinction in- 
creases, we can easily observe reddening effects rising as 
the mean magnitude shifts are more pronounced in bluer 
bands. In addition, we can observe that the mean mag- 
nitude in the z band is almost unaffected by reddening 
in the range considered here. 
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Fig. 1. — SDSS quasar detection sensitivity as a function of reddening. The solid-line histograms show the magnitude distribution of 
SDSS quasars without strong Mgll absorber. The filled histograms show the fraction of objects that remain detectable after reddening 
effects. The vertical lines show the mean magnitude of each sample. 
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Fig. 2. — (Left) Top panel: amount of extinction as a function of E(B — V) for an SMC reddening curve. Middle panel: measured mean 
magnitude shift of the detected objects. We observe that the measured extinction is lower than the initially induced one. Note that the 
mean z band magnitudes are not significantly affected by exti nction. Lower panel: fraction of SDSS quasars still detectable after having 
applied a given extinction. Similarly to the analysis done by Richards ct al. (2003), the SDSS quasar algorithm is able to recover most 
quasars reddened by E(B — B) < 0.1. 

Fig. 3. — (Right) Measured reddening as a function of induced reddening. Even if the measured extinction differs from the induced ones, 
the color changes are not affected by the corresponding selection effects and in the reddening range of interest for the present study, the 
observed color changes of quasars can directly be used in order to quantify the amount of induced reddening by intervening absorbers. 
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Fig. 4. — left panel: redshift distributions of detected Mgll absorbers with Wo > 0.8 A (filled histogram) and their background 
quasars (solid line). Middle panel: rest equivalent width distribution of the Mgll systems detected in SDSS quasar, right panel: observed 
distribution of relative rest equivalent width uncertainty. 



We present the summary of the completeness and red- 
dening results in Figure [2j From top to bottom we show, 
as a function of E(B — V), the input extinction, the mea- 
sured extinction and the fraction of objects that can be 
detected. These results show that the observed extinc- 
tion, defined as the difference between the mean mag- 
nitude of the reddened quasars and those of the initial 
sample, is always smaller than the induced extinction. As 
explained in section [4.1.51 such an effect is expected and 
is due to the shape of the quasar luminosity function. 
Moreover, we can see that the observed mean z-band 
magnitude is not significantly affected by extinction ef- 
fects (for the range of reddening values of interest in this 
analysis). This property allows us to use the z band to 
constrain gravitational lensing effects. The bottom panel 
of the figure allows us to quantify the fraction of objects 
missed as a function of observed E(B — V). By knowing 
the amount of reddening induced by a class of absorber 
systems one can use this curve to infer the fraction of 
objects missed by the SDSS quasar finder. 

In the right panel of Figure [2j we show another 
property related to the selection of reddened quasars: 
the induced and measured colors are the same in 
the range of reddening values of interest in our case, 
i.e. with E(B — V) < 0.5. This indicates that the 
fraction of quasars lost due to color cuts is small. 
This property therefore allows us to measure the 
mean color excess of a sample of quasars and use the 
corresponding color value to infer the induced extinction. 

The properties and selection biases of the SDSS tar- 
get quasar algorithm will be used below to interpret the 
measurements. Our analysis is restricted to the main 
sample of SDSS spectroscopic quasars (with i < 19.1) 
in order to quantify certain selection effects. Including 
fainter quasars is expected to change the above num- 
bers by a modest amount as the slope of quasar number 
counts dist ribution is a slowly va rying function of mag- 
nitude (see iRichards et all (|2005l )). Redshift-dependent 
selection effect s are n ot considered here. As shown by 
IRichards et all (gll), the SDSS quasar selection algo- 
rithm is less efficient in detecting reddened quasars at 
z ~ 2.2. This effect can be neglected in the present anal- 



ysis as only a small fraction of objects are found at such 
a redshift (see Fig. [4J. 



3.3. Absorber systems 
3.3.1. The Mgll absorber sample 

In order to investigate the effects of reddening and 
magnification due to intervening absorbers, we use the 
sample of Mgll system s comp iled with the method pre- 
sented in INestor et ail (|2005f ) and based on SDSS DR4 
data. In this section we briefly summarize the main steps 
involved in the absorption line det ection procedure. Fo r 
more details we refer the reader to INestor et al.l ([2005) . 

All quasar spectra from the SDSS DR4 database 
were analyzed, regardless of QSO magnitude. The 
continuum-normalized SDSS QSO spectra were searched 
for Mgll AA2796, 2803 doublets using an optimal extrac- 
tion method employing a Gaussian line-profile to mea- 
sure each rest equivalent width Wo- All candidates were 
interactively checked for false detections, a satisfactory 
continuum fit, blends with absorption lines from other 
systems, and special cases. The identification of Mg II 
doublets required the detection of the 2796 line and at 
least one additional line, the most convenient being the 
2803 doublet partner. A 5<t significance level was re- 
quired for all A2796 lines, as well as a 3cr significance 
level for the corresponding A2803 lines. Only systems 
0.1c blue- ward of the quasar redshift and red- ward of 
Ly— a emission were accepted. For simplicity, systems 
with separations of less than 500 km/s were considered 
as a single system. The line finder is able to detect Mgll 
absorption lines with a rest equivalent width Wo > 0.3 A. 
However a high completeness is reached only in the range 
Wo > 0.8 A. We will therefore focus our analysis on these 
stronger systems only. In such a range, multiple ab- 
sorbers detected in the same quasar spectrum are found 
only in ~ 15% of the systems. For such cases we consider 
only the stronger system. We note that Mgll absorption 
lines are in general saturated and in these cases no col- 
umn density information can be directly extracted from 
Wq 2796 . The redshift distribution of the systems used 
in this analysis is presented in Figured The filled his- 
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togram shows that of the Mgll systems and the solid line 
their background quasars. The middle panel of Figured 
shows the observed distributions of absorber rest equiva- 
lent width Wq. The intrinsic distri bution, i.e. correcte d 
for incompleteness is presented in iNestor et al.l ([2005) . 
The right panel of the figure shows the distribution of 
rest equivalent width measurement uncertainty. It shows 
the existence of a correlation between line detectability 
and absorber rest equivalent width. As can be seen, the 
relative Wo uncertainty becomes larger for weaker sys- 
tems. Because we require the 2796 A absorption lines 
to be detected at least at the 5-sigma level, this results 
in a selection bias: weaker lines are preferentially found 
in high S/N spectra. In addition, given the fact that all 
SDSS quasars are taken with a similar exposure time, the 
selection criteria introduce a correlation between quasar 
brightness and spectrum S /N. The line finder is therefore 
expected to preferentially find weaker systems in brighter 
quasars. A crucial step in our analysis will be to properly 
take this effect into account when we define a reference 
population of quasars without absorbers. 

3.3.2. The simulated Mgll samples 

In order to test the efficiency of our detection tech- 
nique and to identify biases and systematic effects we 
ran Monte Carlo simulations of the absorber catalog to 
create a sample of simulated Mgll absorption lines re- 
producing the same properties as the real sample (z, Wq, 
A Wo, etc.) and same biases, but with the difference that 
the absorption lines are put in random spectra, i.e. that 
no correlation exists between quasar magnitude and ab- 
sorber properties. Such a population can then be used 
as a control sample. 

In order to create this sample, we proceed as follows: 
for each detected system, a simulated doublet having the 
same redshift, and similar Wo (2796), doublet ratio, and 
Gaussian-fit width was put into many randomly selected 
EDR spectra. Each spectrum containing a simulated 
doublet was then run through the entire non-interactive 
and interactive pipelines, and z and Wo (2796) were mea- 
sured for detected systems. Over 9100 simulated dou- 
blets were put in random quasar spectra. A large frac- 
tion of them appear in regions of spectra with insufficient 
signal to noise ratio to meet the detection threshold for 



the input Wo (2796). The line finder is then able to re- 
cover about 4400 systems. Measurement error of the 
rest equivalent width Wo (2796) can cause systems with 
Wo (2796) < Wiim to scatter into the output catalog as 
well as systems with Wo (2796) > Wu m to scatter out. 
Therefore, the recovered distribution of rest equivalent 
width differs from the input one. In order to create a sim- 
ulated absorption line sample with statistically identical 
distributions of rest equivalent widths, a trial Wo(2796) 
distribution was specified and then adjusted so that the 
simulated output best represents the data. We initially 
chose lines randomly from the input catalog according to 

Wq 

a distribution of the form N(W (2796)) oc exp w «« with 
an initial guess for Wi n , until the number of lines recov- 
ered were equal to that of the actual catalog. Lines with 
input W (2796) < 0.3 A were used, although only lines 
with recovered Wo (2796) > 0.3 A were retained. We de- 
termined a maximum-likelihood value for W ou t using the 
recovered Wo (2796) values. We then corrected our guess 
for Win to minimize \W out - W da ta\- 

This procedure allows us to create a simulated absorp- 
tion line sample having the same properties as the real 
sample described in the previous section. The distribu- 
tions of real and simulated absorber redshifts, rest equiv- 
alent widths and rest equivalent widths uncertainties are 
presented in Fig[5J As can be seen, the statistical prop- 
erties of the simulated samples match very well those of 
the real data. The figure shows that weaker absorption 
lines are preferentially found in high spectrum S /N - and 
therefore brighter quasars - in the same way as the real 
systems are. This bias is systematically introduced by 
the requirements of the line finder algorithm. Below we 
show that the selection of reference quasars described in 
the previous section takes this effect into account. 

3.4. Defining reference quasars 

As our goal is to isolate the magnitude shifts induced 
by intervening absorber systems, the selection of refer- 
ence quasars must (i) not depend on quasar magnitudes 
and (ii) take into account the fact that the absorption 
line detectability depends on the spectrum S/N as seen 
above. 

Here we describe how we define such populations of 
quasars: we first select the quasars with detected ab- 
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sorbers above a given rest equivalent width Wq 
0.8 A. We will compare them to quasars without any de- 
tected absorber above this limit. In the goal of avoiding 
biases that would give rise to magnitude shifts, a number 
of criteria must be applied: the two quasar populations 
must have the same redshift distribution and present the 
same absorption-line detectability. For a given absorber 
with rest equivalent width Wq and redshift z a b Sl we ran- 
domly look for a quasar without any detected absorber 
(above W A2796 = 0.8 A) such that: 

(i) they have similar redshifts, i.e. z Q soref =z QS oabs± 

0.1 ; 

(ii) in addition, we need to check whether this quasar 
would have allowed for the detection of the absorption 
line, i.e. if its S/N at the corresponding wavelengths 
is high enough. Such a requirement translates into 
W m in(QSO re f, z ahs ) < Wo, where W min is the minimum 
equivalent width that can be detected by the line finder 
at the corresponding wavelength. The quasars selected 
in this way are called reference quasars in the following. 
Our selection procedure ensures that if similar absorbers 
population were present in front of the selected reference 
quasars, we would have been able to detect them. 
Since the number of available quasars without detected 
absorbers is larger than the number of quasars with ab- 
sorbers for the range of Wo we are interested in, we repeat 
the previous step Nb times in order to create a series of 
Nb reference samples that will be used for estimating the 
noise associated to a single quasar sample. 

4. ANALYSIS & RESULTS 

In order to isolate the magnitude changes due to the 
presence of the absorber systems we measure, in a given 
band j, the mean magnitude difference between quasars 
with absorbers and reference quasars: 

A {raj) = (TO abs ,j) - • (10) 

As explained above, we can repeat the selection of ref- 
erence quasars Nb times and create a series of reference 
samples in a bootstrap manner. The intrinsic dispersion 
of the series samples can then be used in order to esti- 
mate the errors in the measured magnitude shifts. As 
an example, we show in Figure [6] the results obtained 
by using a sample of quasars with Mgll absorbers hav- 
ing 1.4 < Wo < 2 A and their corresponding reference 
quasars. The left panels show the magnitude distribu- 
tions of quasars with absorbers (solid line histograms) 
and reference quasars (filled histograms). Within the 
Poisson noise, no significant difference can be seen bet- 
ween the two types of distributions. In the right panels, 
we show the mean magnitudes of these two populations. 
For each band, the vertical line shows the value (m a b s ) 
and the histogram shows the distribution of the mean 
magnitudes (m re f) of Nb = 100 reference samples. Such 
an estimator allows to reveal a significant difference be- 
tween the magnitude distributions of the two populations 
and indicates that quasars with absorbers are affected by 
reddening effects. 

Before attempting to detect various signals with the 
data, we make use of the sample of simulated absorbers 
in order to test whether our procedure introduces any 
selection bias. Indeed, as measuring the quantity A (rrij) 
involves the use of various procedures performed over a 
large number of objects (continuum normalization, line 
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Fig. 6. — The left panels show the magnitude distributions 
of quasars with absorbers (solid line histograms) and reference 
quasars (filled histograms). Within the Poisson noise, no signif- 
icant difference can be seen between the two types of distributions. 
In the right panels, we show the mean magnitudes of these two pop- 
ulations. For each band, the vertical line shows the value (m a b s ) 
and the histogram shows the distribution of the mean magnitudes 
( m ref ) °f Nb reference samples. 



detectability, rest equivalent width fitting, etc.), it is im- 
portant to investigate the presence of systematic effects 
and quantify them. 

To check whether the properties of quasars with (de- 
tected) random absorption lines match those of quasars 
without, we proceed as follows: for a given sample of 
simulated absorbers, we define a population of reference 
quasars as defined in section T3.4I We then compare the 
mean magnitude difference between the two populations. 
Since the simulated absorption lines were initially put in 
random spectra, any departure from zero will quantify 
the level of systematics associated to our selection pro- 
cedure. Such a process can therefore test whether the 
samples of quasars with and without absorption lines to 
have the statistical properties. We measure the corre- 
sponding mean magnitudes and the results are summa- 
rized in the left panel of Figure [3 

As can be seen, the mean magnitudes of these two 
samples are similar, but a non-zero magnitude difference 
is found and does not seem to depend on absorber rest 
equivalent width or wavelength. Using the z-band as 
a reference (the least affected by extinction effects, see 
section l4.1.5j) . the mean offset is found to be 

A(m) bms = -0.06 ±0.01 mag. (11) 

Other bands give similar results, ranging from 
-0.05 ± 0.01 to -0.06 ± 0.01. This negative value 
indicates that our global procedure, i.e. line finder 
and/or reference quasar selection, preferentially selects 
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Fig. 7. — Mean magnitude shifts induced by the presence of absorber systems, as a function of rest equivalent width. Left: the results 
obtained using simulated absorption lines in random quasars show the existence of a bias giving rise to an excess brightening of Am ~ 0.06 
(indicated by the gray region). Right: Mgll absorbers detected in SDSS quasar spectra induced significant reddening and extinction effects. 
When the selection bias is taken into account, no significant quasar brightening can be observed. 



absorbers in brighter quasars. This systematic effect 
occurs at the 6% level in terms of quasar magnitude 
or similarly spectrum S/N. It could be corrected by 
increasing the estimates of Wu m by a similar amount. 
Here we prefer to take this effect into account by consid- 
ering the existence of an offset in the mean magnitude 
differences. It should also be pointed out that our 
procedure indicates that the selection of absorption lines 
does not introduce any color-dependent bias. Indeed, 
no significant excess can be detected by comparing the 
colors of quasars with simulated absorption lines to 
those of reference quasars. 

We now apply the same procedure to the sample of real 
Mgll absorbers. The results are presented in the right 
panel of Figure [7l The gray region shows the zero point 
offset, A(m)bias, given by the Montecarlo simulations. 
As we can see, it appears very clearly that the presence of 
a Mgll absorber reddens the light of background quasars 
and the stronger the absorber rest equivalent width, the 
redder the quasar. The direct measurements indicate 
that A(m) < for the strongest systems. However, by 
taking into account the offset due to the selection bias 
A (m) bias, no significant difference can be observed. This 
result contrasts with a number of previous studies claim- 
ing that absorber systems are preferentially found in the 
spectrum of brighter quasars. It suggests that previous 
claims of such correlations may have arisen due to selec- 
tion biases. The current analysis shows that the use of 
Montecarlo simulations in order to test the recovery of 
random absorption lines is an important task to achieve 
in order to quantify the amount of systematic effects. 
More details on the reddening and magnification effects 
are presented below. 

4.1. Dust reddening 

From the magnitude shifts measured in each band we 
measure the corresponding color excess induced by Mgll 
systems. The measured color excess induced by the pres- 
ence of absorbers provides us with constraints on the 
properties of the dust associated with these systems such 



as the size distribution of the dust particles and the re- 
lated column densities. 

It should be noted that the magnitude shifts measured 
in different bands (fig. [7|) are correlated and the errors 
on the colors are in general significantly smaller than 
the errors on the magnitude shifts. We have therefore 
computed the color covariance matrix from bootstrap re- 
sampling and estimated the errors on the measured color 
changes. The results are presented in fig. [8] where we 
show the mean color excess measured for the four color 
combinations u— g, g—r, r—i and i—z, as a function of ab- 
sorber rest equivalent width. As can be seen, strong Mgll 
absorbers induce, on average, reddening values ranging 
from E(g — r) oos ~ 0.02 to ~ 0.3, which correspond to 
similar E(B — V) b s values. 

As the SMC extinction curve is known to describe, 
on average, the reddening properties of Mgll absorber s 
(jKhare et alJ2005tlMinard et al.l20"05l : lYork et alJ2006f >. 
we can write 

(E(B-V))~(N H {z)k{z)) x [£(A B )-£(Ay)](12) 

where Nh{z) is the hydrogen column density associated 
to a Mgll system, k(z) is the dust to gas ratio, and £(A) 
is the wavelength dependence of the extinction curve. 
This relation shows that investigating the color excess 
as a function of redshift allows us to constrain the mean 
product (Nh(z) k(z)). In Figure [8] we have plotted the 
SMC-based reddening values for a series of amplitudes 
normalized by the SMC dust column density 



N d 



N, 



d.SMC 



(io^) cm 



(13) 



where k is the dust-to-gas ratio and Njj is the hydrogen 
column density. Our calculation includes the redshift 
distribution of the corresponding absorbers in each bin. 
First, we can see the SMC-like dust provides a good fit to 
the measured reddening for most of the absorber systems. 
The dotted lines can be used to quantify the amount of 
dust among various strengths of Mgll absorbers. They 
correspond to N^/N^smc — 1,3,10,30 and show that 
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the optical depth for extinction spans a factor ~ 30 for 
Mgll systems in the range Wq = 1 to 6 A. 

If the dust-to-gas ratio of Mgll absorbers is compa- 
rable of that of the SMC, it implies that the strongest 
Mgll absorbers have, on average, hydrogen column den- 
sities of Njj ~ 10 21 cm -2 , i.e. correspond to the regime 
of damped Lyman-a systems. This result is in agreement 
with the hydrogen column densities of Mgll-selected sys- 
tems measured by Rao et al. (2005). 
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Fig. 8. — Color changes of quasars due to the presence of Mgll 
absorbers, as a function of absorber rest equivalent width. The 
solid lines indicate the expected color changes due to an SMC red- 
dening curve with different normalizations. The match between 
the slope of this curve and the data points show that the dust 
related to Mgll systems is similar to SMC dust. 



As the sensitivity of the SDSS filters is significantly 
lower in the u and z bands compared to the <?, r and i 
ones, the g — i color is measured the most accurately and 
will be used below to quantify reddening properties. We 
then convert the color excesses E(g — i) into E(B — V) 
and present the color excesses as a function of absorber 
rest equivalent width. The results are presented in Fig. 
El As can be seen, we find a strong correlation between 
these two parameters. We can also observe that the re- 
lation between observed mean E{B — V) and Wo is well 
described by a simple power-law: 



■;E(B-V).,,,.,; (ir l) ) c 



(14) 



with C = (0.8 ± 0.1) x 10~ 2 , a = 1.88 ± 0.17. As shown 
in section [27X1 for an SMC-type extinction law, the color 
excess E(B — V) is roughly proportionally to the optical 
depth for extinction. Hence we can also write 



(<7y)(Wo) = r Vfi 



Wo 
lA 



(15) 



with ry i0 = (2.5 ± 0.2) x 10~ 2 . 

The relation holds for systems with 1 < Wo < 7 A 
which correspond to velocity dispersions spanning 



100 < Av < 700 km/s and the dust column densities 
change by a factor ~ 30 in this range. The weakest 
and strongest absorbers of this sample are therefore 
significantly different but our results indicate that 
their dust and gas content follows a universal relation, 
suggesting that one phenomenon regulates, on average, 
the amount of one quantity with respect to the other. 
Modeling the gas distribution in the halo of galaxies is 
beyond the scope of this paper but we note that this re- 
lation strongly constrains theoretical models of absorber 
systems, as the dependence relates spatial information, 
Tv, to a quantity that mostly carries velocity space 
information, Wo- We believe that understanding this 
relation will greatly improve our knowledge about the 
nature of these systems. 

If we consider the global distribution of strong Mgll 
absorbers it is interesting to estimate the relative con- 
tribution of optical depth for extinction dC T /dWo, as a 
function of absorber rest equivalent width. Such a quan- 
tity is given by the product: 



ON 



dWo dW 



x (tv)(Wo). 



(16) 



The number of absorbers per unit rest equivalent width 
and redshift is given by Nestor et al. (2005): 



dN/dW6 



A2796 



N* 
W* 



e 



(17) 



with the maximum likelihood values W* — 0.702 ± 
0.017 A and N* = 1.187 ± 0.052. This estimation takes 
into account incompleteness detections and is valid down 
to a rest equivalent width limit of ~ 0.3 A. Using this 
result, we have computed the relative contribution of 
dust dC T /dWo, as a function of absorber rest equiva- 
lent width and present the results in the upper panel of 
Figure [3 The data points show that, for Mgll systems 
with Wo > 1 A, most of the dust contribution originates 
from systems with Wo ~ 1 — 2 A. Therefore most of the 
extinction effects are expected to be associated to such 
systems. It is interesting to mention that an excess of 
Mgll absorbers with Wo — 2 A is reported along gamma 
ray burst lines-of-sight with respect to those of quasars 
(Prochter et al. 2006). 

We can also attempt to estimate the relative dust con- 
tribution originating from all Mgll absorbers, i.e. in- 
cluding also weaker systems. First, we know that the 
functional form used for the incidence of Mgll absorbers 
dN/dWo 2796 holds down to systems as weak as ~ 0.3A. 
In the range 0.3 to 1 A we do not have any detection of 
dust reddening but we can attempt to extrapolate the 
scaling found above (eq. \W\i . By doing so, we find 
that the global dust contribution originating from sys- 
tems weaker than 1 A is smaller than ^20%. The result 
is similar if, instead of extrapolating eq. 1161 we use an 
upper limit from systems with Wo — 1 A, i.e. the weak- 
est systems used in the present analysis. Therefore, the 
vast majority of the dust probed by Mgll absorption lines 
originates from strong systems, with Wo > 1 A. As dust 
originates from stars, this result appears to be in line 
with the association between strong Mgll absorbers and 
~ L+ galaxies, i.e. galaxies from which most of the star 
light originates. 
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Fig. 9. — Left: Mean E(B-V) values observed as a function of absorber rest equivalent width Wo- A power-law fit gives: (E(B — V)) = 
(0.008 ± 0.001) X Wq 88 ^ - 17 . Right, top panel: the relative contribution of mean dust column density, dC T /dWo, as a function of Mgll 
rest equivalent width. The dashed line shows the behavior obtained by extrapolating the observed relation between E(B — V) and Wo- 
Using such an assumption, we show the cumulative distribution of dust column density C(< Wo) in the bottom panel. Most of the dust is 
expected to be carried by Mgll systems with Wo > 1 A. 



4.1.1. Dust redshift evolution 

We now investigate the dust column densities associ- 
ated with Mgll systems as a function of redshift. Using 
an SMC-type extinction curve, we have computed the 
rest frame E(B-V) values as a function of rest equivalent 
width for three redshift bins. The results are presented 
in Figure [TOl and show that the scaling between E(B — V) 
and Wo found in the previous section is consistent with 
being redshift independent but the overall amplitude de- 
creases with redshift indicating that Mgll absorbers are 
less dusty at earlier times. 

In order to quantify the redshift evolution of the dust, 
we measure the average E(B — V) of all absorbers sys- 
tems with Wo > 1 A, in three redshift bins spanning 
the range 0.4 < z < 2.2. In each case we measure the 
observed and rest frame E(B — V) values. The results 
are shown in the right panel of Figure 1101 As we can 
see, the observed reddening effects of Mgll systems are 
on average consistent with being redshift independent. 
Our vision of the distant universe is equally affected by 
Mgll absorbers at z ~ and z ~ 2. In contrast, the rest 
frame E{B — V) values, and therefore the dust column 
densities, show a significant redshift evolution. We find 
that the amount of dust associated with Mgll absorbers 
increases by more than a factor two between z = 2 and 
z = 0.4. A simple power law fit to the data points gives 
E{B-V) reat x{l + z)- 1 - 1±0A . 

It is interesting to see that such a variation is simi- 
lar to that of the cosmic star density fl*. We illustrate 
this agreement in the figur e by showing the integrated 
star formation rate given bv lHernquist fe Springe! |2003) 
where the scaling has been renormalized to the redden- 
ing values. Our results therefore provide a new and in- 
dependent method to probe the star formation rate over 
cosmological time scales. 



If we assume that the redshift and rest equivalent 
widths dependencies can be separated (as suggested by 
the data points), we can write that 

(E(B - V) Test ) (Wo, z) = C X (j^j (1 + zf , (18) 

where a = 1.88 ± 0.17 (see section OJ, = -1.1 ± 0.4 
and C = (0.60 ±0.07) x 10~ 2 . This functional form pro- 
vides us with interesting constraints on the properties of 
cold gas and dust. It will be particularly interesting to in- 
clude such information in theoretical models of absorber 
systems. 

4.1.2. Dust-to-metals ratio 

We can use the above results to estimate the dust-to- 
metals ratio of Mgll absorbers by computing the ratio 

Kdm = (E(B - V))/{N{ZnII)) , (19) 

adopting Zn (an undepleted Fe-peak element) as an in- 
dicator of overall metal content. This ratio provides 
us with an estimate ot the fraction of refractory ele- 
ments depleted onto dust grains, a quantity determined 
by the balance between the formation and destruction 
of dust particules. Using a similar sam ple of Mgll ab- 
sorbers with z ~ 1. iNestor et al. I (|2003f) have measured 
the mean rest equivalent width of Zn//A2026 and found 
Wo = 0.033 ±0.005 A. Assuming that most of the zmc is 
in the form of Znll, we obtain N z „ = (3.63 ± 0.5) x 10 12 
cm~ 2 . Combining this value with our estimate of the 
mean E(B-V) of Mgll systems at z ~ 1 we find 

Kdm = (4.1 ±0.6) x 10" 15 magcm 2 . (20) 

As shown by Wild et al. (2006, see their table 5), the 
dust-to-metals ratios, 72-tjm, of the Milky way, LMC and 
SMC are 4.7, 3.1 and 1.6 x 10~ 15 mag cm 2 , respectively. 
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Fig. 10. — Left: Rest frame E(B-V) induced by Mgll absorbers measured for three redshift bins, and as a function of Mgll rest equivalent 
width. Right: Evolution of the global E(B-V) (for all Mgll systems with Wo > 1 A) as a function of redshift. This trend shows that the 
the mean product (Nh(z) k(z)) increases by a factor ~ 2 between z = 2 and z = 0.5. This is comparable to the evolution of the cosmic 
star density (or integrated star formation rate) illustrated here by the (renormalized) model by Hernquist and Springel (2005). 



Our analysis therefore shows that, on average, Mgll ab- 
sorbers with Wq > 1 A have a dust-to-metals ratio sim- 
ilar to that of our Galaxy. The properties of as a 
function of redshift and Mgll rest equivalent width will 
be presented in a future paper. 

4.1.3. Dust extinction bias in magnitude limited QSO 

samples 

From the observed E(B-V) values, we can now use 
the results of the reddening Montecarlo simulations pre- 
sented in section 13.21 and infer the fraction of quasars 
with absorbers missed in the magnitude limited sam- 
ple. Note that our procedure takes into account quasars 
missed due to both magnitude and color selections. As- 
suming that the distribution of intrinsic reddening values 
is not strongly skewed (which is suggested by the results 
of lEllison et alJ (|2004l ) with radio-selected quasars), we 
estimate the fraction of Mgll absorbers missed due to 
extinction effects and present the results in Fig. [11] The 
fraction of missed systems is below one percent for ab- 
sorbers with Wo — 1 A. However, above this value, in- 
completeness effects become more important. Our anal- 
ysis indicates that the SDSS misses more than 10% of 
Wq ~ 3 A Mgll absorbers and this rises to more than 
50% at Wq ~ 6 A. By extrapolating the trend to stronger 
systems, we see that the SDSS could not detect Mgll 
systems with Wq > 10 A. The fraction of missed quasars 
with absorbers can be parametrized by 



(fm) (Wo) = a 



Wo 

lA 



(21) 



with C m = 0.005 ± 0.001 and 7 = 2.23 ± 0.18. 



The knowledge of f m allows us to estimate the global 
fraction of quasars with strong Mgll absorbers, per unit 



redshift, missed due to the extinction bias. We have 

f°° BN 
f* = l U^^dWo 

r °° dN 

/m(Wo)— dWo 

^0.01 



dW 



(22) 

Therefore, if we consider quasars at z = 1, about one 
percent of them cannot be observed due to extinction 
effects by Mgll systems. As both dN/dWo and the ob- 
served E(B — V) of Mgll systems do not vary strongly 
with redshift (see section |4.1.1|) . this missing fraction is 
roughly proportional to the quasar redshift. Therefore, if 
we consider the main sample of SDSS quasars (for which 
the mean redshift is above 1.5) the incompleteness due 
to the dust hosted by Mgll absorbers is about 2%. 

4.1.4. Recovering the unbiased dN/dWo distribution 

As the fraction of missed quasars with absorbers de- 
pends on Wo, the observed distribution of absorber rest 
equivalent widths differs from the intrinsic one and we 
have 



dN 
dWo~ 



( dN 



\dW 



1 



l-/m(Wo) 



(23) 



As we have used the functional form of dN/dWo to es- 
timate the missing fraction f m , an exact estimation of 
these quantities should be made with an iterative pro- 
cess using eq. [23] However, given the small values of f m 
and the exponential cutoff of dN/dWo, such an approach 
will only bring negligible corrections to f m . Our previous 
estimations are therefore robust. 

The values of f m shown in Figure [Tl] indicate that 
the shape of the observed distribution of dN/dWo is ex- 
pected to be significantly affected at the high-end values 
and to present a deficit. 
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brightening effect, which differs from a number of 
previous studies claiming detections of magnification 
effects due to metal absorbers. It is interesting to point 
out that without the quantification of systematic effects 
done with the Montecarlo simulations, i.e. without 
the calibration of the zero point, the raw measurement 
with the real data would have led to the detection 
of a negative magnitude shift which could have been 
interpreted as magnification effects. 
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Fig. 11. — Fraction of quasars with Mgll systems that is ob- 
served due to reddening and extinction effects. For the strongest 
systems we can see that almost half of the objects are missed in 
this magnitude limited sample. 



4.1.5. Gravitational lensing 

In order to detect the potential effects of gravitational 
lensing, we now estimate the magnitude shifts induced by 
the presence of strong Mgll absorbers, after correcting 
for dust extinction. 

In section 13.21 we have shown that, on average, the 
induced and measured color excess are equivalent (in 
the range of values of interest in our analysis). How- 
ever, given the shape of the QSO magnitude distribution 
and the selection criteria of the SDSS quasar target algo- 
rithm, a difference exists between induced and measured 
extinction. Having quantified these effects in section I3~2l 
we can now use the observed E(B — V) values and ap- 
ply the appropriate observed extinction correction to the 
quasar magnitudes. The corresponding results are shown 
in Fig. [T2] where the horizontal line represents the zero 
point obtained with the Monte Carlo simulations for sim- 
ulated absorption lines. By comparing the extinction- 
corrected magnitude shifts to this reference level, we do 
not find any significant magnification signal (for which 
Am would be negative). 

If we consider all strong Mgll absorbers with Wo > 
l A, we find the mean magnitude shift in the z band to 
be: 

(Ait) = 0.028 ±0.015 (24) 
(Ag) = 0.015 ±0.014 
(Ar) = 0.015 ±0.013 
(Ai) = 0.015 ±0.013 
(Az) = 0.010 ±0.013 

which was compared to the zero-point and taking into 
account its uncertainty. These mean magnitude shifts 
are consistent with zero ((An) is positive at the 2-a 
level only). Our analysis does not show any detectable 
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Fig. 12. — Mean QSO magnitude shift due the presence of a 
strong Mgll absorber, after correcting for dust extinction, as a 
function of absorber rest equivalent width Wo ■ The horizontal line 
represents the zero point calibrated with Monte Carlo simulations 
using simulated absorption lines in real SDSS quasar spectra. The 
colors denote the filter names (see Fig. 0. Our analysis does not 
show any significant magnification effect, which would correspond 
to a negative (Am) 



We can use the above results to infer an upper limit on 
magnification effects. If we consider the g, r and i bands 
the 3-cr upper limit on magnification is 



(Am) 



lensing 



-0.024. 



(25) 



As explained in section 14.1.51 and detailed in Menard 
(2005), given the shape of the quasar magnitude distri- 
bution, this corresponds to an upper limit for the mean 
magnification: 

(M)<1-10, (26) 

meaning that the excess magnification is smaller than 
about ten percent. In order to convert this value into a 
more interesting quantity, we can assume the mass dis- 
tribution of the galaxies responsible for the absorption to 
follow that of a singular isothermal sphere. The galac- 
tic halos therefore have the following density and surface 
density: 



p ex 



and 



2Gr 



(27) 
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and the corresponding magnification effects of a point 
source can then be computed in a simple way: 

_ / 2/y if y < 1 



/' = 



1 + 1/y if y > 1 



(28) 



where y is the impact parameter normalized to the Ein- 
stein radius of the lens: 

m 

where a v is the velocity dispersion and -Dd,s,ds & re the an- 
gular diameter distances from the observer to the lens, 
to the source, and from the lens to the source. Using the 
above upper limit on the magnification, and the average 
impact parameter for Mgll systems with Wo > 1 A re- 
cently measured by Zibetti et al. (2006): (b) ~ 40 kpc, 
we find 

cr„ < 269 kms" 1 , (30) 

for absorbers at z = 0.5 and quasars at z — 1. This result 
shows that, on average, Mgll absorbers are not associ- 
ated with massive elliptical galaxies and is in agreement 
with previous estimates based on different methods. Pre- 
vious works have shown that the average Mgll absorb- 
ing galaxy has a luminos ity about ~ 0.8 L* ()Zibetti et al.l 
120071 iSteidel et al J 119971) . The expected velocity disper- 
sion of the gravitational potential of such galaxies is ex- 
pected to be ~ 150 km s _1 . Therefore the sensitivity 
of the present measurement needs to be improved only 
by a factor of a few in order to detect the magnification 
signal originating from the galaxies associated with Mgll 
absorbers. The present analysis is based only on quasars 
with i < 19.1. In the future, it might be of interest to 
quantify the selection effects of fainter SDSS quasars as 
it would double the sample size. In a longer term, being 
able to measure this magnification signal as a function 
of absorber rest equivalent width and redshift will pro- 
vide us with important constraints on the nature and 
environment of metal absorbers. 

5. SUMMARY 

We have used a sample of almost 7000 strong Mgll ab- 
sorbers with 0.4 < z < 2.2 detected in the SDSS DR4 
dataset to investigate the gravitational lcnsing and dust 
reddening effects they induce on background quasars. 
We have attempted to make significant improvements in 
making such measurements, by using a sample of Mgll 
absorbers an order of magnitude larger than earlier ana- 
lyzes and by carefully quantifying a number of systematic 
effects previously neglected. 

In order to do so we have restricted our analysis to 
a well defined sample of ~ 30,000 SDSS quasars with 
i < 19.1 and considered strong Mgll absorbers with a 
rest equivalent width Wo > 1 A. We have quantified the 
efficiency of the absorption-line finder algorithm with 
Monte Carlo simulations and the sensitivity of the SDSS 
quasar target algorithm to detect reddened quasars. 
We then measured the statistical magnitude changes 
induced by the presence of Mgll absorbers. 

Our main results are as follows: 

(i) Strong Mgll absorbers significantly redden the light 
of their background quasars. We confirm previous re- 
sults showing that, on average, the dust particle size 



distribution of these systems is similar to that of the 
SMC, i.e. without the presence of a reddening excess at 
0.2/im. Moreover, we find that the observable reddening 
effects follow the simple relation: (E(B — V) bs) (Wo) = 
C (W ) Q , with C = (0.8±0.1) x 10~ 2 , a = 1.88±0.17, for 
absorbers with 1 < Wo < 6 A. In this range, the mean 
dust column density increases by a factor ~ 30. Since 
most Mgll lines with Wo > 1 A are saturated, this scal- 
ing suggests that, on average, the dust column density 
is proportional to the velocity dispersion of the gas Such 
a relation provides us with an important constraint for 
theoretical models of absorber systems. 

(ii) The analysis of reddening effects as a function of 
has shown that the rest-frame E(B — V) values follows 
(E{B~V)) oc (\ + z)-f\ with (3 = 1.1±0.4, which implies 
that the amount of dust in Mgll systems increases by a 
factor ~ 2 between z = 2 and z — 0.5. Such a trend is 
similar to the evolution of the cosmic star density. Our 
results therefore provide a new and independent method 
to probe the star formation rate through the production 
of dust over cosmological time scales. 

The observed E(B — V) values induced by Mgll appear 
to be consistent with no redshift evolution down to z — 2. 
This observation results from less dust at higher redshift 
being canceled out by higher SMC extinction at high 
redshift (which probe bluer rest-frame wavelengths). We 
have also showed that, considering all Mgll absorbers, 
most of the dust is carried by systems with Wo ~ 1 — 2 A. 

(iii) We have estimated the dust-to-metals ratio of 
Mgll absorbers by measuring the quantity TvLdm = 
{E(B -V))/{N Z „) = (4.1 ±0.6) x 10- 15 magcm 2 which 
is similar to that of the Milky way. 

(iv) We have quantified the fraction of absorbers 
and/or quasars missed due to extinction effects induced 
by Mgll systems. While less than 2% of absorbers with 
W ~ 1 A are missed, this fraction increases up to ~ 30% 
for the strongest systems of our sample, with Wo ~ 6 A. 
This effect therefore affects the shape of the distribution 
of absorber rest equivalent widths and we provide a cor- 
rection factor for it. We also show that extinction effects 
due to Mgll systems decrease the number of observable 
quasars by less than 2%. 

(v) Regarding gravitational lensing, in contrast to 
previous studies we do not find an excess of absorbers in 
brighter quasars and we obtain only an upper limit on 
gravitational magnification: /j,(Wq 796 > lA) < 1.10. We 
have shown that testing the global selection procedure 
(for detecting absorption lines and defining reference 
quasars) is crucial in order to make such a measurement. 
Our result provides an upper limit for the velocity 
dispersion of Mgll absorbing galaxies: a v < 269 km 
s" 1 . As it has been shown that the mean luminosity of 
Mgll absorbing galaxies is about ~ 0.8 L+ (|Zibetti et al.l 
2007), it implies that improving the sensitivity of our 
measurements by a factor of a few might lead to the 
detection of the magnification effects induced by Mgll 
absorbing galaxies and therefore a characterization of 
their mass. 
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6. FUTURE DIRECTIONS 

Understanding the nature of the structures probed by 
absorption lines is an important task to achieved. Among 
metal lines, the Mgll doublet, AA2796,2803A is one of 
the most commonly detected features in optical quasar 
spectra but little is known about the properties and ori- 
gin of Mgll absorber systems. Attempts to find correla- 
tions between Wo and galaxy parameters have turned out 
to be rather difficult. Recently, using H ST images of 37 
galaxi es giving rise to Mgll absorption, iKacprzak et al.1 
(|2007h showed that most galactic parameters do not cor- 
relate with absorptioin strength. 

In this paper we have shown the existence of a well- 
defined scaling relation between the mean E(B-V) in- 
duced by Mgll absorbers and their rest equivalent width 
Wo, and presented its evolution as a function of redshift. 
Since most Mgll lines with Wo > 1 A are saturated, this 
scaling translates into a relation between dust column 



density and gas velocity dispersion, with Av spanning 
the range 100 to about 600 km/s. Theoretical models 
are needed in order to gain insight into the origin of this 
relation. We believe that being able to understand and 
reproduce this scaling will shed light on the nature and 
properties of the structures traced by low-ionization ab- 
sorption lines. 
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